Shape of the Polymerization Rate in the Prion Equation. 
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Abstract 

We consider a polymerization (fragmentation) model with size-dependent parameters involved 
in prion proliferation. Using power laws for the different rates of this model, we recover the shape 
of the polymerization rate using experimental data. The technique used is inspired from [30] where 
the fragmentation dependency on prion strains is investigated. Our improvement is to use power 
laws for the rates whereas |30j used a constant polymerization coefficient and linear fragmentation. 
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1 Introduction 

Polymerization of prion proteins is a central event in the mechanism of prion diseases. The cells 
produce naturally the normal form of this protein (PrPc) but there exists a pathogenic form (PrPsc) 
present in infected cells as polymers. These polymers have the ability to transconform and attach the 
normal form by a polymerization process not yet very well understood. Mathematical modelling is a 
promising tool to investigate this mechanism provided we can find accurate parameters. 
We consider the following model of prion proliferation (see [U [2j [10] [12] ) 



0_ 

at 



dVjt) 
dt 

u(x, t) 

u(0,t) 
u(x, 0) 



A - V(t) 



8 + 



t(x)u(x, t) dx 



d 

-V(t)—(r(x)u(x,t)) - [/j,o + P{x)]u(x,t) + 2 j P(y)K(x } y)u(y,t)dy, (i) 
0. 

U Q (x). 



This equation is the continuous version of the discrete model of [17] (see [3] for the derivation). 
System (JTJ) is a quite general model able to reproduce biological experimental results (see [15]). In this 
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equation, V(t) is the quantity of PrPc which is produced by the cell with a rate A, and degraded with 
the rate 5. The quantity of aggregates of PrPsc with size x at time t is denoted by u(x,t). Polymers 
of PrPsc can attach monomers (PrPc) with rate r and can split themselves in two smaller polymers 
with rate (3. To fit with experimental observations [29], the polymerization and fragmentation rates 
should depend on the polymer size x (see [HE]). When fragmentation occurs, a polymer of size y 
produces two pieces of size x and y — x with the probability n(x,y). The natural assumptions (see 
[20l EU E5 CO]) on the kernel k are then 

Supp k = {(x,y), < x < y} 

since fragmentation produces smaller aggregates than the initial one, 

Vy > 0, J K.(x,y)dx = 1 («;(■, y) is a probability measure) (2) 

and 

Vx < y, n(x,y) = n(y - x,y) (3) 

because of indistinguishability of the two pieces of sizes x and y — x after fragmentation. A consequence 
of the two assumptions ([2]) and (|3|) is that the mean size of the fragments obtained from a polymer of 
size y is y/2 

Vy, J xn(x,y)dx = |. (4) 

In the discrete model of [W\ \T7\ I18j . a critical size n of polymers is considered. Under this size, 
the polymers are unstable so the proteins that compose them go back to their normal form PrPc. 
This effect is taken into account in the continuous model of [10t [TT] with the critical size xq. As 
a consequence, there is one more integral term of apparition of monomers in the equation on V. 
Nevertheless, even if n > in the discrete model, the value of xq can be taken equal to in the 
continuous one as justified in [3]. Then the fact that small polymers are unstable, and so do not exist, 
is taken into account by the bounday condition u(0,t) = 0. The parameter fj,Q is a death term which 
will be discussed later. 

In the article [30], model ([I]) with the specific parameters of P [TUl El HI \17\ [18] is used. It assumes 
that the polymerization rate does not depend on polymer size (r(x) = To) and the fragmentation is 
described by (3(x) = /3qx and n{x,y) = \< x <y\- In this case, equation (fTJ) can be reduced to a non- 
linear system of ODEs (see [Hllinilll] or even [T6J [TTl [18] ) . This system is used in [30] to investigate the 
parameters involved in the so-called "prion strain phenomenon" . This phenomenon is the ability of the 
same encoded protein to encipher a multitude of phenotypic states. Prion strains are characterized 
by their incubation period (time elapsed between experimental inoculation of the PrPsc infectious 
agent and clinical onset of the disease) and their lesion profiles (patterns produced by the death of 
neurones) in brain [9]. A reason for the strain diversity can be that there exist various conformation 
states of PrPsc, namely various misfoldings of the prion protein, which lead to various biological 
and biochemical properties \12\ I27j. These different conformations could be characterizable by their 
different stability against denaturation. Indeed, recent studies on prion disease have confirmed that 
the incubation time is related not only to the quantity of pathogenic protein inoculated (J uo(x) dx) 
and the level of expression of PrPc (A), but also to the resistance of prion strains against denaturation 
|14j in terms of the concentration of guanide hydrochloride (Gdn-HCl) required to denaturate 50% of 
the disease-causing protein. Other studies have highlighted a strong relationship between the stability 
of the prion protein against denaturation and neuropathological lesion profiles |13|. 128] , A natural 
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question is then to know how the conformation states of the PrPsc protein affect the parameters 
of model ([I]) and more precisely which of these parameters is the most strain-dependent. One of 
the conclusions of [30] is that the key parameter describing strain-dependent replication kinetics is 
the fibril breakage rate (5q. Using the model, some coefficients already defined in [T7J are calculated 
from experimental data and the dependency between them indicates that the fragmentation rate /3q 
is the most adapted parameter to be strain-dependent. But, in this framework, the fragmentation 
dependency is not sufficient to fit precisely the data and the polymerization rate To is also assumed to 
change with the strain. In particular, to keep a one-parameter family of models, the polymerization 
rate is assumed to depend on the fragmentation one through a power relation 

30 s.t. r = (5) 

However, polymerization and fragmentation are two different biochemical processes which have no 
reason to be linked together. 

In this paper we prove that we can get rid of such a relation by considering a larger class of parameters 
for model ([I]) . We take advantage of a general framework developped in [H El H2] , where coefficients 
can be more general than in [6j [TUl [HI [El [T71 [181 ED] , and we consider power laws as in [71 [121 [19] . 
Thanks to this extension, we find a shape for the polymerization rate allowing us to fit experimental 
data with only the fragmentation rate dependent on strain. The experimental data we use are the 
stability against denaturation G and the growth rate of the quantity of PrPsc (r) which is linked to 
the incubation time. These data are strain-dependent and we can find in [30J estimated empirical 
values for different strains. In Section [2] we present the class of coefficients we consider and give the 
necessary assumptions for mathematical investigations. Most of these assumptions are already made 
in |17[ [30] and, thanks to them, we are able to link r and G to the mathematical model. In Section [3] 
we give the dependency of these parameters on polymerization and fragmentation rates. Then we give 
a relation between them using the only fragmentation dependency on strain and are able to fit data 
with an accurate shape of the polymerization rate. 



2 Assumptions and Method 

To avoid the use of a power relation like ([5]) between the polymerization and the fragmentation, we 
allow parameters r and f3 to depend on the size x with power laws : 

37, v and /3o, To s.t., Vx > 0, [3(x) = /3o^ 7 and t(x) = tqx u . (6) 
Concerning the fragmentation kernel k, we assume that there exists a measure kq on [0, 1] such that 

/ x \ 

Vx,y>0, K(x,y) = -K - . (7) 



y 

Then the natural assumptions on k become 



kq(z) dz = 1 and kq(z) = kq(1 — z). (8) 







We assume that the exponents v, 7, the kernel kq and the death rate /j,q are strain-independent unlike 
the "intensities" /3q and tq. First we look which of /3q and tq is the most characteristic of the strain. 
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Then we find a value of v allowing to fit experimental data as well as [30j but without assumption ([5]). 
We will see in Section [3] that the method naturally gives the parameter v because it turns out that 7 
and Ko do not appear in the final argument. 

The method proposed in [TTJ [30] is to investigate the dynamics of the system (pQ) linearized at the 
disease-free equilibrium. Then we can use the eigenelements as suggested in [TJ [2]. To make this 
approximation valid, we make some classical assumptions. First we assume that, in the biological 
experiments (namely before the animal death), the amount of PrPc in cells remains close to the sane 
value V = j (obtained by setting u = in the first equation of ([I])). Quantity V being assumed to be 
constant, we consider the eigenelements r and U of the equation for polymers associated with V 



rU(x) = -V^ {t(x)U(x)) - [ Mo + P{x)]U(x) + 2 / P(y) K (x, y) U{y) dy, 



U(0) =0, J U{x)dx = 1. 



(9) 



The principal eigenvalue r represents the growth rate of the solutions and the eigenvector U is the 
corresponding size distribution of polymers. To consider that the size distribution is close to this 
profile along the experiments, we assume that the initial value uq{x) is already close to U. It is a 
reasonable assumption if contaminations are made from other infected individuals. Then the solution 
of (HJ) is reduced to 

V(t)~V and u(x,t) ~ [ [ u (y)dy] U(x)e rt . (10) 



To solve our problem, we use the idea that the eigenelements r and U change with the prion strains 
(see [2j [15] for instance) . Because present experiments do not allow one to measure the size distribution 
U, we link it to the stability against denaturation G by the relation 

3a > 0, b s.t. G = ax + b (11) 

where x = J xU{x) dx is the mean size of the polymers since J U(x) dx = 1. Such an affine dependency, 
suggested in [30], is based on the idea that the longer x is, the larger G will be. In [T3J [30] we can 
find estimated empirical values of G (as a concentration, in Mol, of Gnd-HCL required to denaturate 
50% of PrPsc) and of the growth rate r (computed from measures of incubation times thanks to the 
assumption of an exponential growth homogeneous in the brain). These values are reported in Table [U 



Prion strain 


139A 


ME7 


BSE 


Sc237 


RML 


vCJD 


Chandler Scrapie 


301 V 


T(day~ L ) 


0.05 


0.024 


0.015 


0.11 


0.18 


0.07 


0.17 


0.07 


G(M) 


2 


2.9 


2.8 


1.6 


1.7 


1.85 


2.2 


2.2 



Table 1: Estimated empirical values of the growth rate and the stability against denaturation for different 
prion strains (data from (T3l [30] ) . 



In the next section, we find a relation between r and IA by investigating their dependency on 
parameters. Then we show that this relation allows to fit the data of Table [U with an accurate choice 
of the polymerization exponent v. 
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3 Results 

Proving the existence of solutions to the eigenvalue problem ([9]) is not easy and the theory was 
first developed in |25} I19j. In [5] we could treat singular polymerization and fragmentation rates as 
investigated here. It is proved in [5j [19] that solutions exist under the condition 

7 + l-i^>0. (12) 

Moreover, following the idea of [8], we are able to give the dependency of r and IA on parameters 
V, no, and tq for given 7, v and kq (unknown but strain- independent). Define r\ and IA\ by the 
eigenequation given for /3q = Vtq = 1 and no = by 

nWiOz) = - - s?Wi(s) + 2 /~yW-Wy) — . (13) 

Then rx, U\ depend only on 7, v, kq and we have the following lemma which we use later on. 
Lemma 1 The eigenelements r and U solution of © satisfy 



r=(Vr ) T ^Po T ^r 1 -Ho and W( I ) = (^) 1 ^((^) 1 ^ 1 )' 
Proof. The proof is based on a homogeneity argument. First we write ([9]) as 



(14) 



U(x) = ~(x"U(x))-£-x*U(x) + 2p- r y^KJ-)u{y)^- (15) 
ox Vtq Vt J x \yJ y 



Vr 

and we set a = y^r and R = T y^ ■ The rescaled function v defined by U(x) = a k v(a k x) satisfies the 
equation 

a k Rv{x) = -a k ^- v ^{x u v{x))-a k+1 -^x^v{x) + 2a k+l ~ k ^ [°° y^K (- )v(y) — . (16) 
ox J x \yJ y 

We choose k = 1+ }^_ v (what is possible since (fT2|) is satisfied) and obtain 

ex. i+7— " Rv{x) = —-^-(x v v(x) s j — x r v(x) + 2 f y^Ko (— }v(y) — . (17) 
ox j x \y/ y 

Comparing f)13|) and (|17j) . the uniqueness of eigenelements ensures that 

1-1/ 

R = a 1 +t-"n and = Wx(x) (18) 

because f v(x)dx = J U(x)dx, and the lemma is proved. □ 

We are now able to give the dependency of r and s on parameters (3q, tq and no- The mean length 

/ xu(x , i) dx 

of the aggregates at time t is defined as the quotient ^ — - — ^— - — . So, using the approximation (iTOl) . 

J tt(x, t) ax 

we obtain that this quantity is constant in time and equal to 

xU(x)dx. (19) 
Then, thanks to Lemma [U we conclude to the following corollary. 
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Corollary 1 Under assumptions ©-([8]) and (|12p . there are two constants r\ and x± (depending only 
on 7, v and kq) such that 

r = [Vt ) t +^ fo^^rx - /i and x = 0^) x±. (20) 



Proof. Using Lemma [H we immediately have x = ^ — ° ^ 1+7 " x\ with x\ = J xlA\ 



(x) dx. 



□ 



We have obtained the dependency of r and x on the a priori strain-dependent parameters p\) and 
tq. Using assumption (jlip . we deduce the dependency of the stability against denaturation G 

G = a(^) TT ^x- 1 + b. (21) 

In paper [30| . the exponent 7 of the fragmentation rate is equal to 1. Here we do not impose any 
value for 7 but, nevertheless, we consider that it is positive. Indeed the larger the polymers are, the 
easier they break themselves. Consequently a strain-dependency of tq cannot lead to a decreasing 
correspondance between r and G, as observed in Figure [Q since we have assumed that 7 + 1 — v and 
a are positive. Thus it indicates that the key strain-dependent parameter should be /3o- With this 
assumption, we obtain an implicit expression of G as a function of r 

G = A(r + ^ )^+b (22) 
where A and b are constants (namely strain- independent). This model is used to fit the data. 



5 r 
4.5 - 

4 - 
3.5 - 

3 - 
2.5 

2 - 
1.5 - 

1 - 

0.5 - 

- 




Q(r)=A.(r+n„; 



~0 



0.1 

Rate of growth (r) 





1 

Ar*- 1 +b 


A(r+fj, )—+b 


V 


-0.482 


0.316 







0.023 


A 


0.083 


0.01 


b 


1.54 


1.69 


H 2 


0.7 


0.72 



Figure 1: Experimental values of G and r for different Table 2: Best fitting parameters obtained with the 
strains are plotted (circles). Then these data arc fitted two different models. The correlation coefficient R 2 is 
using the reduced model (f22|) with /^o = and uq free, also reported in the last row. 



First we assume as in [30] that the death rate /io is equal to zero. Then we obtain a model already 
considered by [30] (Table 3, page 5). In our framework, we can also let fiQ be free and fit the data with 
the model (|22p . As we can see in Table [2] the accuracy of the fitting is a little bit better with this 
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generalisation than in the case = 0. But the most significant difference between the two fittings is 
the shape of r : it is a decreasing function of x in the first case [y < 0) and it is an increasing function 
in the second case (y > 0). 



4 Conclusion and Perspectives 

Thanks to a larger class of parameters than [30] , we are able to fit the experimental data by considering 
only the fragmentation rate to be dependent on strains. We even obtain a more general model and 
so a more accurate fitting. But, because the parameters profiles are not a priori known, we cannot 
find x from experimental data and have to postulate a particular dependency (llip whereas this was 
a consequence of the model in |30j . To validate it, new data are necessary, for instance values of x 
obtained by an experimental size distribution. 

A consequence of the method is the prescription of the shape of the polymerization rate from 
experimental data, and also the death rate fiQ. Nevertheless, we notice a significant difference between 
the two recovered exponents v obtained by assuming that jjlq = or not. Particularly one is negative 
and the other one positive (see Table [2]). The idea to compare parameters from different strains 
has been deeply exploited here. But the fact that there are only few different strains in comparison 
to the number of fitting parameters is a weakness for precise fitting and precise determination of 
r. Some parameters such as /j,q or b should be previously determined by another method and other 
data. Furthermore the parameters of model ([TJ) could have a size-dependency more general than power 
laws. A complete inverse problem is to recover all the size-dependent parameters from experimental 
data without the restriction to a class of shape. This problem is very complicated even though first 
elementary models are treated in [H [261 E3]- Moreover it requires more information such as the size 
distribution U and not only the mean size x. 
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